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ABSTRACT 


The  equations  implicit  in  Noel's  molecular -field  theory  of 
ferrimagnetism  are  set  out  in  considerable  detail.  The  general 
equations  which  pertain  to  magnetic  structures  with  crystallo¬ 
graphic  long-range  order  are  derived,  and  the  specific  pair  of 
transcendental  equations  which  must  be  satisfied  in  the  case  of 
two  sets  of  antiparallel  magnetic  spins  is  written  down  in  a  form 
which  is  amenable  to  solution  by  automatic  computer  methods. 
The  dominant  interaction  is  expressed  in  terms  of  the  magnetic 
transition  temperature,  while  the  two  smaller  interactions  are 
given  as  fractions  of  the  dominant  one.  A  computation  of  the 
shapes  of  the  curves  of  the  magnetization  as  a  function  of  tem¬ 
perature  shows  that,  while  magnetic  and  neutron  diffraction  meas¬ 
urements  can  determine  quite  precisely  the  value  of  a  linear  func¬ 
tion  of  these  fractions,  the  actual  values  of  the  individual  fractions 
are  less  precisely  determined.  It  follows  that  any  experimental 
deviation  from  the  calculated  curves  is  due  to  an  inadequacy  of 
the  theory  rather  than  to  incorrect  values  of  the  parameters. 
Finally,  methods  of  attack  for  computation  of  more  complex  cases 
are  discussed. 
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CALCULATION  OF  MAGNETIZATION  VS  TEMPERATURE  CURVES 
FOR  FERRIMAGNETIC  MATERIALS 


INTRODUCTION 

The  N£el  theory  of  ferrimagnetism  (1,2)  is  recognized  universally  as  providing  a 
good  basis  for  understanding  the  magnetic  properties  of  a  wide  variety  of  substances 
which  contain  magnetic  ions  on  two  or  more  crystallographically  distinct  sets  of  sites  in 
their  crystal  structure.  The  theory  basically  assumes  that  it  is  possible  to  describe  the 
total  magnetic  energy  of  the  crystal  by  an  expression  of  the  type 

E  =  -L  Ju  s,  -s, .  (1) 

where  the  summation  is  taken  over  all  pairs  of  ions  i  and  j ,  and  S,  and  Sj  are  the 
angular  momentum  vectors  of  the  ith  and  jth  atoms  respectively.  The  quantity  JM  is 
the  exchange  integral,  which  is  a  measure  of  the  strength  of  the  interaction  tending  to 
align  the  angular  momenta  of  the  two  ions.  If  J,  j  is  positive,  the  energy  will  be  lowest 
if  S,  and  s;  are  parallel  to  one  another,  and,  conversely,  a  negative  value  of  J, }  favors 
antiparallel  alignment  of  S,  and  S, .  Under  the  assumption  of  Eq.  (1)  the  contribution  to 
the  total  energy  due  to  the  ith  ion  alone  has  the  same  dependence  on  the  angular  orienta¬ 
tion  of  S,  as  it  would  have  if  the  ion  were  a  magnetic  dipole  acted  upon  by  a  magnetic  field. 
Equation  (1)  can  then  be  written  in  the  form 

'  ‘ ?  s‘ '(£ "“)  <»• 
and  the  summation  over  j  is  termed  the  "molecular  field"  acting  on  the  ith  ion. 

In  most  crystals  we  may  assume  that  it  is  possible  to  divide  each  set  of  crystallo¬ 
graphic  sites  into  a  finite  number  of  subsets  in  each  of  which  the  time-averaged  spin 
directions  are  all  parallel,  and  we  also  assume  that  the  individual  contributions  to  the 
molecular  field  are  proportional  to  the  expectation  values  of  gSj ,  where  g  is  the  gyro- 
magnetic  ratio.  Now 


is  defined  as  M, ,  the  magnetization  of  the  set  of  ions  over  which  the  sum  has  been  taken. 
Therefore,  we  may  express  the  molecular  field  acting  on  the  ith  ion,  or  on  any  other  ion 
in  the  same  subset  in  the  form 


(3) 


The  coefficients  a. }  are  the  "molecular  field  coefficients"  for  the  interaction  ij .  The 
values  of  these  coefficients  determine  the  shapes  of  the  curves  of  the  magnetizations  |M,  j 
as  a  function  of  temperature.  These  individual  magnetizations  are  often  called  the  "sub¬ 
lattice  magnetizations,"  although  the  sets  of  crystallographic  sites  are  often  not  lattices. 
These  curves  can  be  measured  by  means  of  neutron  diffraction,  and  sometimes  by  other 
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methods,  and,  at  least  in  principle,  the  values  of  the  molecular  field  coefficients,  and 
therefore  an  approximation  to  the  values  of  the  exchange  integrals,  can  be  inferred  from 
the  shapes  of  the  measured  curves. 

The  purpose  of  this  report  is  to  derive  the  mathematical  expressions  for  these 
curves  in  a  form  amenable  to  computation,  to  set  down  the  specific  equations  which 
apply  to  the  simple  case  of  two  sets  of  ions  with  their  magnetic  moments  aligned  anti¬ 
parallel  to  one  another,  to  describe  the  results  of  calculations  made  using  these  equa¬ 
tions,  and  to  indicate  the  possible  procedures  in  extending  the  calculations  to  more  com¬ 
plex  situations.  The  appendix  describes  in  detail  a  program  written  for  the  NAREC 
computer  to  do  the  calculations  in  the  two-set  case. 


GENERAL  EXPRESSIONS  FOR  THE  MAGNETIZATION 
VS  TEMPERATURE  CURVES 


It  can  be  shown  (3)  that  an  ion  in  a  magnetic  field  H  will  have  an  expected  value  of 
its  magnetization  given  by  M  =  g  S Bs(x),  where  g  is  the  gyromagnetic  ratio,  is  the 
Bohr  magneton,  and  Bs(x)  is  the  Brillocin  function, 


Bs(x) 


2S  +  1 
2S 


ctnh 


"( 2S  +  l)x  " 
2S 


2S 


ctnh 


(4) 


where  x  =  gS|H|  ^/kT  and  S  denotes  the  magnitude  of  the  vector  S.  The  molecular  field 
approach  assumes  that  the  same  expression  can  be  used  if  H,  instead  of  being  an  applied 
field,  is  the  fictitious  "molecular  field." 


We  shall  resolve  the  system  of  vector  equations  (3)  into  sets  of  three  scalar  equa¬ 
tions  involving  components  along  three  suitable,  orthogonal  directions,  which  will  be 
designated  by  a  running  subscript  k.  Thus, 


H 


ik 


(5) 


If  we  denote  the  direction  cosines  of  the  molecular  field  Hi  by  Cik,  and  remember 
that  Mi  is  parallel  to  H; ,  Eqs.  (5)  become 


(6) 


But  In, |  =  kTx,/g jSj/Xg,  and  |M,I  =  nigjSiMBBsifXj) ,  where  ni  is  the  number  of  ions  per 
formula  unit  in  the  ith  set.  Therefore, 


kT x i  C i k  =  Mbj  8 jS 4  Zj  cij  .  Cjk 


gj  sr 


(7) 


If  there  are  n  sets  of  ion  sites,  there  will  be  3n  equations  like  (7).  A  spontaneous  mag¬ 
netic  moment  can  exist  if  there  exists  a  nonzero  vector  x  with  components  x,  clk  which 
is  a  solution  of  Eqs.  (7),  subject,  of  course,  to  the  restriction  that 


1  . 


We  shall  prove  later  that,  in  certain  special  cases  at  least,  a  nonzero  solution  is  possi¬ 
ble  below  a  certain  limiting  temperature,  which  can  be  identified  with  the  Curie  temper¬ 
ature,  Tc.  Further,  we  can  infer  from  the  laws  of  thermodynamics  that  if  at  least  one  of 
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the  interactions  is  nonzero  for  each  value  of  i  (or  j)  there  must  be  a  solution  for 
very  small  values  of  T,  and  this  solution  must  approach  x,  =  <»  for  all  values  of  i  as  T 
approaches  zero. 


THE  CASE  OF  TWO  SETS  OF  IONS 

The  special  case  in  which  the  magnetic  structure  consists  of  only  two  sets  of  mag¬ 
netic  ions  oriented  antiparallel  to  one  another,  with  a  net  moment  resulting  from  either 
unequal  numbers  of  ions  or  unequal  magnetic  moments  per  ion  on  the  two  sets  of  sites, 
is  of  considerable  practical  importance,  as  well  as  mathematically  tractable.  In  this 
case  we  can  take  one  of  the  component  directions  parallel  to  the  moments,  so  that  Cu  =  1, 


Cj,  =  -1,  ci2  =  Ci3  =  0.  Then  Eqs.  (7)  reduce  to  two  only: 

kTxj  =  /j.g  SjSj  [f>ta  j  jgjSjBjjtXj )  -  njai  282B2Bs2(x2)]  (8^) 

-kTXj  =  mbJ  82B2  [niaij8iSiBsl(x,)  -  n2aJ2g2S2BS2(x2)]  ,  (8b) 

The  equations  can  be  rewritten  in  the  slightly  more  convenient  form 

ntan*isiBsi(*i)  “  njau82SjBsj(x2)  '  (kT/MBJ8isi>  x,  =  0  (9a) 

ia  1  2®  1B  1BS if X  1 )  ^  ^ 2^2 2®* 2B 2BS 2^ x j)  “  (kT//^Bg2S2)  xa  “  (9b) 


We  are  particularly  interested  in  the  case  where  all  the  atJ  are  negative  and  a12 
has  a  considerably  larger  magnitude  than  04,  or  aJ2.  Under  these  conditions,  Eq.  (9a) 
defines  xt  as  an  implicit  function  of  x,  which  is  continuous  and  positive  for  all  positive 
values  of  x2.  Its  derivative  is 


dxl  -  n2al282B2B  S2(x2)  (JO) 

dx2  (kT/ gjSj)  -  nja^gjSjB’sjCx,) 

where  B's(x)  denotes  (d/dx)Bs(x).  For  aI2  and  a,2  negative  this  expression  is  posi¬ 
tive  and  finite  for  all  values  of  x2  provided  either  T  is  greater  than  absolute  zero  or 
an  is  not  equal  to  zero.  At  x  =  0,  B's(x)  =  (S+  l)/3S,  so  that 


dx,  _  _ -n2al  282^S2  +  1V3 _ 

dx2  * 2*o  (kT/MB  g^j)  -  njijjgjCS,  +  l)/3 

This  expression  approaches  zero  as  T  becomes  very  large,  and  has  the  value 

^  2^  1 2®  2^  ^  2 

n,aiigi(S,  +  1) 


(11) 


at  T  =  0. 

By  a  similar  argument,  Eq.  (9b)  defines  x2  as  an  implicit  function  of  x,  for  all  the 
same  conditions.  It  follows  that,  provided  (“,2/an)  >  (“22/ai2),  there  must  exist  a  non¬ 
zero  solution  to  Eqs.  (9)  if  T  is  not  too  large.  It  also  follows  that,  for  T  sufficiently 
large,  the  only  possible  solution  is  x,  =  x2  =  0.  We  may  define  the  Curie  temperature, 
Tc,  as  the  least  upper  bound  of  values  of  T  for  which  a  nonzero  solution  exists.  A  neces¬ 
sary  condition  for  this  temperature  is  that  the  two  curves  defined  by  Eqs.  (9)  should  be 
tangent  at  the  origin,  which  implies  that 
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kT 

"1*1111*8,(8,+ D/3  -  — | 


n2ajl*lB,Sl(S2  + 


kT. 


nla12®2Bl^2^1  +  ^^3  n2a2  2&  2  ^2^2+  l)/3  “  r 


=  0 


(12) 


The  Curie  temperature  is  usually  an  easily  measured  quantity,  so  it  is  convenient  to 
express  the  other  interactions  in  terms  of  it.  To  this  end,  let  us  define  /3,  =  au/ai2, 

/32  =  a j 3/ a j  2 ;  then  we  can  rewrite  Eq.  (12) 


nAkiJsi(si+  D/3  -  kTc  /mbJ  a  j  j  njgjgjS^Sjt  l)/3 

n.gjRjSjtS,  +  l)/3  n2£2g2S  S2(S2  +  D/3  -  kTc/.uB  a,2 


0  . 


(13) 


Then 


kV^B2al2  =  |(n,/3ig,JS1(S1+ 1)  +  n/jg/SjfSj+l)  -{[n1/3lg,2S1(SJ+  1) 

-  nj/3jgjJ  Sj(Sj  +  1)J  2  +  4n,nJgI2gJJS1S2(S1  +  1)(S2+ 1)JI/2  j  , 

where  the  minus  sign  is  chosen  in  the  quadratic  formula  since  Tc  must  be  positive  but 
a12  is  negative.  If  we  denote  the  right-hand  side  of  Eq.  (14)  by  -G,  we  get 


1/a  I, 


and  Eqs.  (9)  become 


~n|/3jgj  Sj  9gj(Xj)  +  n 2g jg2S jS2Bg 2( x 2 )  GTxj/Tc  -  0 
n  1^18 jS jSjB S|(x j )  -  n2/32g2  S2  BS2(x2)  -  GTx2/Tc  =  0  . 


(15) 


(16a) 

(16b) 


These  equations  can  be  solved  by  means  of  numerical  techniques  in  a  fairly  straight¬ 
forward  manner.  First  we  choose  a  pair  of  values  for  x,  and  x2  which  are  fairly  cer¬ 
tain  to  be  larger  than  the  final  solution.  Using  this  x2  we  solve  Eq.  (16a)  for  x,,  and 
using  x,  we  solve  equation  (16b)  for  x2.  We  compare  the  new  values  of  x,  and  x.  with 
the  previous  ones,  and  if  both  changes  are  less  than  a  predetermined  tolerance  we  nave 
found  the  solution.  If  not,  we  find  the  lines  tangent  to  the  two  curves  at  the  two  new  points 
we  have  found,  solve  the  equations  of  these  lines  for  a  new  xt  and  x2,  and  repeat  the 
process,  which  will  converge  in  most  cases  after  a  few  iterations. 


APPLICATION  TO  LITHIUM  FERRITE 

Rado  and  Folen  (4)  have  treated  special  cases  of  molecular  field  calculations  in  two 
set  systems  in  which  both  sets  of  ion  sites  are  occupied  by  the  ion  Fe*  +  t ,  in  which  case 
Sj  =  S2  =  s/2,  and  g,  =  g2  =  2.  With  these  restrictions,  Eqs.  (16)  reduce  to 


-nl^iBj/2^ Xj)  +  n2®s/ 2( *j)  -  GTXj/Tc  -  0 


(17a) 
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(17b) 


where 


G  -  ~  3q  |ni/3,  +  nj/3j  -  [fn,/3,  -  nj^j)  +  i n a ]  }•  (^®) 

Rado  and  Folen  used  a  manual  technique  to  calculate  the  magnetization  vs  tempera¬ 
ture  curve  for  lithium  ferrite,  LiFe,Oa.  This  substance  crystallizes  in  the  spinel  struc¬ 
ture,  which  is  based  on  a  framework  of  oxygens  in  an  arrangement  which  is  a  small  dis¬ 
tortion  of  cubic  close  packing.  The  cations  are  distributed  among  the  interstices,  some 
in  "tetrahedral"  or  "A"  sites  with  fourfold  coordination  and  some  in  "octahedral"  or  "B" 
sites  with  sixfold  coordination.  In  lithium  ferrite  two  of  the  five  iron  ions  occupy  A  sites. 
The  substance  has  two  arrangements  for  occupancy  of  the  B  sites  by  the  remaining  ions, 
depending  on  the  previous  thermal  history  of  the  sample.  If  the  sample  is  rapidly  cooled 
from  a  high  temperature,  the  remaining  three  iron  ions  and  the  lithium  ion  occupy  the  B 
sites  at  random.  On  the  other  hand,  if  the  sample  is  slowly  cooled,  the  lithium  ions 
occupy  a  particular  subset  of  B  sites  (5)  in  a  manner  which  possesses  long-range  order. 

The  numbers  n,  and  n,  do  not  depend,  at  least  to  a  first  approximation,  on  the  degree 
of  the  B  site  ordering  but  depend  only  on  the  total  proportion  of  magnetic  ions  in  the  two 
types  of  sites. 

Rado  and  Folen  were  able  to  find  values  for  the  ratios  /?,  and  /3a  which  gave  excel¬ 
lent  agreement  between  the  observed  and  calculated  curves  of  net  magnetization  vs  tem¬ 
perature.  Because  the  net  magnetization  is,  however,  a  function  of  the  difference  between 
the  magnetizations  of  the  individual  sites,  there  remains  the  possibility  either  that  the 
solution  is  not  unique  or  that  the  success  of  the  molecular  field  method  in  this  case  was 
largely  fortuitous  and  the  magnetizations  of  the  individual  sites  are  not  adequately 
accounted  for.  That  this  might  be  the  case  is  suggested  by  some  recent  nuclear  magnetic 
resonance  work  (6)  in  which  the  magnetizations  of  both  sites  appeared  to  be  higher  than 
predicted  at  temperatures  near  the  Curie  temperature.  In  an  effort  to  determine  the 
sensitivity  of  the  magnetization  vs  temperature  curves  to  changes  in  the  values  of  the 
molecular  field  parameters,  we  have  calculated  the  curves  for  the  parameters  given  in 
Rado  and  Folen’s  paper,  as  well  as  for  some  other  sets  of  parameters  chosen  so  as  to 
give  the  same  value  of  the  net  magnetization  at  the  temperature  0.75  Tc.  The  programs 
used  for  these  calculations  are  described  in  the  appendix.  The  procedure  was  to  calculate 
the  curves  of  the  two  individual  site  magnetizations  and  the  net  magnetization  using  the 
Rado  and  Folen  parameters,  /3,  =  0.54,  =  0.22.  Then  a  new  value  of  /3,  was  specified, 

and  the  computer,  by  standard  numerical  methods  determined  the  value  of  which  gave 
the  same  net  magnetization  at  T  =  0.75  Tc .  Finally  the  computer  calculated  new  curves 
of  magnetization  vs  temperature.  Curves  were  calculated  for  /3,  =  +0.20  and  for  /3,  =  -0.45. 
The  latter  curves  were  obtained  as  a  result  of  an  inadvertent  error  in  coding  the  problem 
for  the  computer,  but  the  results  were  highly  instructive.  Although  this  corresponds  to  a 
positive  value  for  aM,  in  which  case  a  solution  has  not  been  previously  shown  to  exist, 
the  solution  actually  does  exist,  and  in  fact  the  curve  differs  from  the  original  by  only  a 
few  percent.  Figure  1  shows  the  net  magnetization  as  a  function  of  temperature  for  the 
three  pairs  of  values  =  0.54,  /3,  =  0.22;  *  0.20,  /3,  *  -0.0009;  and  *  -0.45,  03  = 

-0.42.  Figure  2  shows  these  values  of  /*,  plotted  against  /3, ,  and  shows  that  there  is  a 
very  nearly  linear  relationship  between  these  parameters  under  the  assumed  system  of 
constraints.  The  similarity  of  the  curves  of  Fig.  1  is  evidence  that  the  parameters 
and  are  very  strongly  correlated,  and  that  measurement  of  the  magnetization  curves 
is  not  a  very  sensitive  measure  of  the  values  of  these  molecular  field  parameters. 

Table  1  shows  the  A-site  magnetizations  alone,  for  the  two  extreme  pairs  of  param¬ 
eters.  These  curves  are  evidently  also  insensitive  to  the  molecular  field  parameters,  so 
that  if  a  current  neutron  diffraction  investigation  confirms  the  nuclear  magnetic  reso¬ 
nance  curves,  it  will  be  evidence  that  the  agreement  of  the  net  magnetization  data  is 
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Fig.  2  -  Plot  of  the  ratio  PJ  against  the  ratio 
/3j  for  a  constant  value  of  net  magnetization 
at  T  =  0.75  Tc  .  Magnetic  measurements  can 
be  expected  to  locate  the  position  of  this  line 
quite  precisely,  but  the  location  of  the  best 
point  on  the  line  is  very  uncertainly  deter¬ 
mined. 


Fig.  1  -  Net  magnetization  vs  temperature 
curves  for  lithium  ferrite  for  various  pos¬ 
sible  values  of  the  ratios  of  the  molecular 
field  coefficients 


largely  fortuitous,  and  that  the  molecular  field  approximation  is  inadequate  in  this  degree 
of  precision. 


PROCEDURE  IN  CASES  OF  MORE  THAN  TWO  ION  SETS 

In  cases  in  which  the  magnetic  structure  includes  more  than  two  sets  of  magnetic 
ions,  the  possible  structure  types  are  much  more  numerous,  and  the  mathematical  treat¬ 
ment  becomes  correspondingly  more  complex.  In  some  instances  one  of  the  interactions 
dominates  all  the  others,  and  the  equations  can  be  solved  in  a  manner  similar  to  the  one 
used  for  the  simple  two-set  case.  In  order  to  eliminate  the  dominant  interaction  param¬ 
eter  in  terms  of  the  Curie  temperature,  it  is  necessary  to  solve  an  equation  similar  to 
Eq.  (12)  but  of  degree  n,  where  n  is  the  number  of  independent  equations.  This  is  likely 
to  be  difficult  except  in  special  cases,  but  it  does  have  the  advantage  of  including  the 
reduced  temperature,  T/Tc,  as  a  parameter  in  the  equations.  Thus  the  scale  of  the  inter¬ 
action  strengths  is  fixed  by  the  relatively  precise  measurement  of  the  Curie  temperature, 
and  the  shapes  of  the  curves  need  only  be  used  to  determine  ratios. 
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Table  1 

Values  of  the  Predicted  A-Site  Magnetization  as  a  Function 
of  Temperature  for  Various  Values  of  the  Molecular  Field 
Coefficients 


T/Tc 

ma/m0 

/3,  =  0.54, 

/3j  =  0.22 

/3,  =  0.20, 

/8j  =  -0.0009 

■ mm 

— 

1.000 

1.000 

Hgns 

Sill/' 

BBiMK 

0.997 

BBS 

mBm-* 

0.990 

0.3125 

Ir  I . 

0.977 

0.3750 

0.960 

0.958 

0.4375 

■SSI,-, 

0.937 

0.P34 

0.5000 

■H  1... 

0.907 

0.905 

0.5625 

0.873 

0.872 

0.870 

0,6250 

0.830 

0.829 

,  ■  ■ 

0.6875 

0.777 

0.776 

1 

0.7500 

0.712 

0.712 

0.712 

0.8125 

0.632 

0.632 

iHwH 

0.8750 

0.528 

0.529 

0.9375 

0.382 

0.383 

0.385 

In  more  complex  cases,  particularly  where  several  of  the  interactions  have  similar 
magnitudes,  so  as  to  produce  structures  containing  canted  spins,  there  is  no  alternative 
to  a  direct  attempt  to  find  solutions  to  Eqs.  (7).  The  starting  point  is  the  structure  at 
absolute  zero,  as  determined  from  magnetic  and  neutron  diffraction  measurements.  At 
absolute  zero  the  arguments  of  all  Brillouin  functions  become  infinite,  with  the  values  of 
the  functions  approaching  unity.  The  value  of  a  given  Brillouin  function  at  a  lower  tem¬ 
perature  can  be  used  as  a  first  approximation  at  a  higher  temperature,  while  better 
approximations  can  be  found  by  numerical  techniques.  In  this  method  the  predicted  Curie 
temperature  must  be  found  by  a  numerical  determination  of  the  highest  temperature  at 
which  a  nonzero  solution  exists.  It  is  a  serious  disadvantage  of  this  procedure  that  the 
independent  variable  is  absolute  temperature  rather  than  the  reduced  temperature,  T/Tc. 
Thus,  in  order  to  fit  experimental  curves  it  is  necessary  to  calculate  the  entire  curve 
before  adjusting  the  constants  to  correct  the  scale  in  such  a  way  as  to  cause  the  solutions 
to  vanish  at  the  correct  Curie  temperature. 

It  is  probable  that  each  individual  problem  will  have  to  be  considered  on  its  own 
merits.  Many  of  the  more  complex  systems  are  three-set  systems  with  antiparallel 
alignments,  such  as  the  inverted  ferrites  of  iron,  cobalt,  and  nickel,  and  gadolinium-iron 
garnet.  These  would  involve  the  solution  of  a  cubic  equation,  which  can  be  handled  by 
analytic  methods,  although  it  is  not  easy  to  see  in  advance  what  the  best  procedure  is 
going  to  be. 
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Appendix  A 

PROGRAM  TO  CALCULATE  MAGNETIZATION  CURVES 


The  program  used  to  calculate  the  curves  of  magnetization  vs  temperature  for 
lithium  ferrite  is  on  NAREC  tape  4658.  This  program  can  be  used  for  magnetization 
calculation  in  any  compound  in  which  there  are  two  sets  of  magnetic  ions,  and  all  mag¬ 
netic  ions  have  the  same  magnetic  moment.  With  some  simple  changes  the  routine  could 
be  used  where  the  magnetic  moments  on  the  two  sites  are  different.  The  program  con¬ 
tains  within  it  a  subroutine  which  calculates  the  value  of  the  Brillouin  function,  Bs(x), 
and  its  derivative,  B's(x) ,  for  values  of  S  from  1/2  to  9/2  and  for  values  of  x  up  to  1000. 
It  also  contains  a  subroutine  which  solves  the  equation  pBs(y)  +  qBs(x)  -  ry  =  0  for  y, 
given  a  value  of  x.  Tape  4658  occupies  NAREC  locations  3caO-3d5f,  and  it  requires  the 
fixed-point  subroutine  tape  4400-3d70  in  locations  3d70-3fff.  Locations  3d60-3d6f  are 
used  as  working  locations.  We  shall  describe  the  individual  parts  of  the  program  in  detail. 


EVALUATION  OF  THE  BRILLOUIN  FUNCTION 

The  Brillouin  function,  Bs(x) ,  is  given  in  Eq.  (4)  in  a  form  which  is  commonly  writ¬ 
ten  because  of  compactness,  rather  than  for  convenience  of  computation.  For  computation 
purposes  it  may  be  rewritten  in  the  form 

(2S  +  1)  [l  +  e  (JS*,),/S]  (1-  e-*/s)  -  (1  +  e  */s)  [l  -  e  (2s+1)*/s] 

B  (x)  =  - - - r -  , 

2S  [l  -  e'<2S‘,>‘/S]  (i_e- ■«/•) 

which  is  less  compact,  but  which  requires  only  the  evaluation  of  two  exponentials  in  addi¬ 
tion  to  simple  arithmetic  operations,  Similarly,  the  derivative  B's(x)  takes  the  form 

-  (2S  +  1)  2  e(2SM>’/s  (1-  e'*78)2  +  e-*/s  [l  -  e‘  <  2St  1  >  */ s]  2 

B*s(x)  =  - - - r-j - j - , 

SJ  [l  -  e-(2S*, >„/*]’  (l  .  e-x/S)2 

which  is  a  simple  function  of  the  same  two  exponentials.  At  x  =  0  both  of  these  formulas 
produce  indeterminate  forms,  but  a  straightforward,  although  tedious,  application  of 
L'Hospital's  rule  gives  the  results  Bs(0)  =  o  and  B's(0)  =  (S+  i)/3S. 

In  order  to  keep  all  quantities  within  the  range  of  values  permitted  in  fixed  point 
arithmetic  on  the  NAREC,  the  argument  of  the  Brillouin  function  is  scaled  by  1000  (deci¬ 
mal).  The  quantity  2S  must  be  placed  in  the  left  address  position  of  the  first  subroutine 
entry  word.  Since  9/2  is  the  largest  acceptable  value  of  S,  it  is  immaterial  whether  dec¬ 
imal  or  hexadecimal  notation  is  being  considered  by  the  programmer.  The  subroutine 
entry  form  is 


s)  ( 2S)<XU>  woU 
s  +  1)  s  +  1  A  L03d20(3d42)  . 

Control  will  be  returned  to  LOs  +  2.  If  the  entry  is  L03d20,  Bs(x)  will  be  in  the  A  register 
on  return.  If  the  entry  is  L03d42,  B„(x)  will  be  in  the  A  register  and  B's(x)  will  be  in 
the  U  register  on  return. 
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SOLUTION  OF  INDIVIDUAL  EQUATIONS 

The  equation  PBs(y)  +  qBs(x)  -  ry  =  0  defines  y  as  an  implicit  function  of  x  for  all 
nonnegative  values  of  x.  A  routine  which  solves  for  y  as  a  function  of  x  obviously  will 
provide  solutions  of  Eqs.  (17a)  and  (17b)  individually,  given  in  each  case  an  initial  value 
for  the  independent  variable.  We  use  here  the  Newton-Raphson  iterative  method,  whereby 
the  second  guess  y,  is  obtained  from  the  first  guess  y,  by 

pBs(y,)  +  qBs(x)  -  ry, 


A  third  guess  is  obtained  from  the  second,  and  so  forth,  until  the  equation  is  satisfied 
within  a  specified  tolerance.  When  the  proper  value  of  y  has  been  found,  the  derivative 
is  computed  from  the  relation 


dy  qB,s(y> 
dx  '  r  -  pB's<x) 


The  subroutine  requires  five  parameters: 

P0  =  P  x  io  - J , 

P,  =  q«  io1, 

P2  =  r  x  10- 1  , 

P3  =  2S  in  the  left  address  position  , 

P4  =  x  x  IO’ 3  . 

These  are  stored  in  consecutive  locations,  the  location  of  P0  being  supplied  to  the  sub¬ 
routine  in  the  entry.  The  entry  form  is 

s)  (address  of  PQ)<XU>  sU 

s  +  1)  s  +  1  A  L03d00  . 

On  completion  the  control  is  returned  to  LOs  +  2  with  10" 3  y  in  the  A  register  and  10’ 3 
dy/dx  in  the  U  register. 


SIMULTANEOUS  SOLUTION  OF  MOLECULAR  FIELD  EQUATIONS 

The  main  subroutine  makes  use  of  the  two  previous  subroutines  to  solve  simulta¬ 
neously  the  two  equations 

-n,/3,Bs(x)  +  n2Bs(y)  -  GTx/Tc  =  0 
niBs(x)  “  *  GTy/Tc  =  0  , 

where 

f  i  ,/J) 

g  =  -  |n A  +  "A  _  [("A  "  n A>  +  4ninj]  j  (J+  i)/6j . 

Starting  from  an  initial  guess  (x0,y0)  the  point  (x,,y0)  is  found  by  solving  the  first 
equation,  and  the  coefficients  of  the  line 


s 

•f 
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K  ~  xl  dx 

y  -  y0  =  dy  <xi>y«> 

are  determined.  Similarly,  the  point  (x0.y,)  is  found  by  solving  the  second  equation, 
and  the  coefficients  of  the  line 


are  determined.  Finally  the  intersection  of  these  two  lines  is  taken  as  the  next  guess, 
and  the  process  is  repeated  until  the  differences  in  the  x  and  >  coordinates  of  two  sue 
cessive  guesses  are  less  than  the  specified  tolerances. 

The  subroutine  requires  six  parameters: 


P0  =  -Pi  .  with  /St  <  1  , 

pi  =  •  with  <  l  , 

P  j  =  n ,  *  10 ' 1  , 

P3  =  n3  x  10- 1  , 

P«  =  2S  in  the  left  address  position  , 

P,  =  T/Te. 

The  entry  form  is 

*1  (address  of  P0)<XU>  sU 
*  +  1 )  s  +  1  A  L03ca0  . 

The  routine  returns  to  the  main  program  at  LOs  +  2  with  Bs(y)  in  the  A  re  gister  and 
Bs(x)  in  the  U  register. 
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